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Abstract: We present a semi-recursive method for calculating the rational parts 
of one-loop gravity amplitudes which utilises axial gauge diagrams to determine 
the non-factorising pieces of the amplitude. This method is used to compute the 
amplitudes M^-'°°p(1-, 2+, 3+, 4+, 5+) and M1-i°°p(1-, 2+, 3+, 4+, 5+, 6+). 



Keywords: [Models of Quantum Gravity, NLO Computation^ 



Contents 

|l|. Introduction |l| 

|2|. Recursion Q 

1^. Axial gauge diagrammatics |^ 

||. The graviton scattering amplitude Afi-'°°P(l-,2+,3+,4+, 5+) || 

|5|. Conclusions and remarks [14 

|A] . Six-point single-minus amplitude [14 

[B|. Graviton scattering amplitudes [16 

[C|. String-based rules calculation of Afi-'°°P(1 , 2+, 3+, 4+, 5+) [17 

|C.1| Summary of the string-based rules for gravity amplitudes [18 

Application to M^-i°°p(1-, 2+, 3+, 4+, 5+) |G 



1. Introduction 

On-shell recursive techniques have proven very successful in the computation of scat- 
tering amplitudes in gauge theories and in theories of gravity [1-3]- The recursive 
techniques for tree scattering amplitudes make use of both the rationality of the 
amplitudes and their complex factorisation properties. Specifically, in a theory with 
massless states, if we use a spinor helicity representation for the polarisation vectors 
it is possible to write the amplitude entirely in terms of spinorial variables A(A^, A^) 
where the massless momentum of the i^^ particle is A^A^ = (cr^)aQ-fcf 

Within this formalism it is possible to probe the analytic structure of the am- 
plitude by choosing a pair a, b of external momenta and shifting these according 
to 

A"^A'^-;^A^ A^— >A'' + 2A" (1.1) 

where we suppress the spinor indices. The analytic behaviour of the shifted amplitude 
A{z) can then be studied. 

^We use the usual spinor products (ab) = e"'^A°A^, [ab] = e^^^AJiA^, which satisfy (a 6) [6a] = 
{ka + fcb)^ = Saf), chains of spinor products such as \a\b\c) = [a 5] (6c) and [a|Pe/ - |^) = [a|e|6) + 
[a|/|6) H etc and t^hc = (^a + fcfc + fcc)^- 
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If A{z) 



1. is a rational function, 

2. has simple poles at points Zi, and 

3. vanishes as z — > oo 

then applying Cauchy's theorem to A{z)/z with a contour at infinity yields 

A(0) = -5^Res(^,z,y (1.2) 



z 

Zi 



This technique has proven very effective in computing tree amplitudes and has been 
extended from the purely gluonic case to a variety of other applications including 
that of gravity [2,3]. Alternate shifts exist [4] which can be used to re-derive the 
CSW formulation for Yang-Mills [5] and gravity [6] . 

The result (|1.2|) holds even if condition ^ above is relaxed to poles of finite 
order; however this condition allows us to use the factorisation theorems to determine 
the residues in terms of lower point amplitudes. At tree level the factorisation is 
relatively simple: amplitudes must factorise on multi-particle and collinear poles. 
For a partition of the external momenta {Sl, Sr) with at least two momenta on 
either side, and defining = J2ieSL n-point tree amplitude A^^"^^ factorises 

as K becomes on shell as 



A';:i{k, G SL,K^)^AT.W^{i-K)-'^,k e Sr) 



1.3) 



where a denotes the internal state of the intermediate particle and r is the length 
of Sl. Consequently, simple poles in the shifted amplitude A{z) occur at values of 
z where K'^{z) = 0. Only those i^"'s containing precisely one of ka or kh will be z 
dependent. When the corresponding K'^{z) vanishes the residue will be the product 
of the tree amplitudes defined at z = Zi. Thus we can express the n-point tree 
amplitude in terms of lower point amplitudes, 

ATiO) = E<H-T(^.);^<-r/i(^0> (1-4) 

where the summation over i is only over factorisations where the a and b legs are on 
opposite sides of the pole. 

Beyond tree level there are three potential barriers to using recursion. Firstly the 
amplitudes, in general, contain non-rational functions such as logarithms and dilog- 
arithms; secondly, the amplitudes may contain higher-order poles for complex mo- 
menta and, finally, the amplitudes may not vanish asymptotically with z. Nonethe- 
less a variety of techniques based upon recursion and unitarity have been developed. 
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A one-loop amplitude for massless particles may be expressed as 



^i-ioop^ J2 cdl + R + 0{e) (1.5) 

n=2,3,4;i 

where the scalar integral functions are the various scalar box, triangle and bubble 
functions. The function R contains the remaining rational terms. The one-loop am- 
plitude can then be specified by computing the coefficients, Cj, and the purely rational 
term R. The Cj are rational coefficients which can be computed by various applica- 
tions of the four- dimensional unitarity technique [7-9] or indeed recursively [10]. 

There are a variety of strategies for evaluating the rational terms. They may 
be evaluated using dimensional unitarity, by recursion or by specialised Feynman 
diagram techniques [11-24]. In general, the rational term R does not simply satisfy 
the previously-stated requisites for recursion ^ and |^. If the amplitude has only 
simple poles but does not vanish as z — > oo then it can be possible to formulate 
recursion by the use of an auxiliary recursion relation [25] . However there are rational 
amplitudes for which one cannot find a shift which only generates simple poles such 
as the one-minus amplitude y4^"'°°P(l^, 2+, ■ ■ ■ , n+). These amplitudes vanish at tree 
level and consequently are purely rational at one-loop. A shift on these amplitudes 
yields double and single poles 

'^~(7^+(7^ + ---=(7^(' + a^"-"^) + ---) ^'-^^ 

The double pole is not in itself a a barrier to using recursion with the double pole 
contributing 

• i \ ^ 1 

^z{z-Ziy' 7 zf 

However to obtain the full residue in a recursive construction one must have specific 
formulae for this double pole and for the coincident single pole, or the 'pole under 
the double pole'. 

In ref. [26] the form of the pole in Yang-Mills was postulated to be 



+ E s{b,, k-, 62) j 



where the 'soft' factors are 

Sia, s\ b) = ^'^f S(a, 6) = (1.9) 

(a s) {s b) [a s\ [s b\ 

With this ansatz recursion correctly reproduces the known one-minus one-loop am- 
plitudes. In ref. [27] it was shown that the consistency requirements for recursion in 
QCD are sufficient to determine these soft factors. 
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The above postulate, or variations thereof, however does not work for gravity 
amphtudes [28]. In this article we will demonstrate how to apply recursion tech- 
niques in gravity scattering amplitudes by determining the 'pole under the pole' 
using an axial gauge formalism. By only keeping the pole terms it is relatively 
simple to extract these from the diagrammatic approach. We demonstrate this 
by calculating the previously-unknown amplitudes M-'^"'°°p(1~, 2+, 3+, 4+, 5+) and 
M(l~, 2+, 3+, 4+, 5+, 6+). We assume that the shifted amplitudes have vanishing 
behaviour as z — )■ oo. The expressions we derive have the correct symmetries and 
soft limits, providing strong evidence for the validity of this assumption. Further, 
we compare the result numerically with a completely independent computation of 
M"°°P(l-,2+,3+,4+,5+) from 'string-based rules' for gravity [29-32]. 



2. Recursion 



The factorisation of one-loop massless amplitudes is described in ref. 



^i-ioop ^£^0 



+ ■ ■ ■ , ki+r-1, j^A*^^^_^_i(^(—K) ^, ki+r, ■ ■ ■ , ki-i)Fn[K^; ki, . . . , kn) 

(2.1) 

where the one-loop 'factorisation function' F„ is helicity-independent. Naively this 
only contains single poles, however for complex momenta there are double poles. 
These can be interpreted as due to the three-point all-plus (or all-minus) one-loop 
amplitude also containing a pole 



^^-'°°P(i^+,a+,6+) = -^yi-i°°P(i^+,a+,6+) 



where, for pure Yang-Mills, 

y'-'°°p(x+, a+, b+) 



487r2 



[Ka] [ab] [bK] 



(2.2) 



(2.3) 



To see this explicitly, let us consider the five-point Yang-Mills amplitude 
^^-•°°P(l-,2+,3+,4+,5+) [34]: 



A^-'°°P(l-,2+,3+,4+,5+) = 



487r2 (34)^ 



[2 5]^ ^(11)^^ [15] (3 5) (13)-^[3 2](12) 



[12] [51] (12) (2 3) (4 5)' (15) (5 4) (3 2) 



(2.4) 
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If we carry out a complex shift on legs 1 and 5, 



X^^X^ + zX\ X^-^X^-zX\ (2.5) 

then (4 5) — > (4 5) + z (4 1) which vanishes at 2; = —(4 5) / (4 1) and the amplitude 
has a double pole at this point. ^ 

A recursive approach suggests drawing the diagrams shown in figure |I]. The 
third of these involves the one-loop vertex V^~^°°'^{K~^ ,4:~^ ,5^). Computing with this 
does correctly generate the double pole in the amplitude [26,28], however it needs 
augmentation to give an expression with the correct single pole. By trial and error, 
adding the second term in (|1.8|) gives the correct single pole and completes the 
computation of the amplitude. 



rT)^LV4+ 2+^L)^T: 






Figure 1: Diagrams contributing to the recursive construction of A{1~, 2~^, 3"*", 4"'", 5^) 
with legs 1 and 5 shifted in the manner of (2.5). The diagram (c) contains the one- 
loop vertex V^^"'°°P(_ft'"'", 4"*", S"*") that contributes the double-pole. 

When calculating the gravity amplitude mI~^°°^{1~, 2^^, 3^, 4"'', 5"'') we must con- 
sider the same class of diagrams as in figure |I| together with permutations over the 
external legs. For gravity the vertex 

V'-'°"nK\a\b+) = -^^([iTa] [ab] [bK]f (2.6) 

can be used to generate a double pole term which has the correct soft and coUinear 
limit, but attempts [28] to implement a universal correction for the single pole anal- 
ogous to that of (|1.8| ) have failed. 

We find that the resolution is to replace the factorisation term of figure |l^ with 
a tree insertion diagram of the form shown in figure ^ and compute this using axial 
gauge diagrammatics. In section |^ we present the axial gauge rules, in section || the 
computation of the five point one-minus gravity amplitude and in appendix ^ the 
result of the computation of the six point one-minus amplitude. 



^The term which gives rise to the double pole [4 5] / (4 5)^ is the one- loop splitting function [7] 
which only gives a linear coUinear pole for real momenta. 
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3+ ^ 




2+ - 




i- ' 


i-5 



^5+ 

Figure 2: The form of the tree insertion that augments the recursion in order to 
construct the double pole and its underlying single pole. 



3. Axial gauge diagrammatics 

We use axial gauge diagrammatic methods to determine the singular structure nec- 
essary to augment the recursion. We identify and compute the singularities arising 
when we shift a negative-helicity leg a and a positive-helicity leg h as in ( |1 . 1|) . These 
singularities arise from propagators involving just external momenta and from the 
loop momentum integration. 

Following ref. [35] we use a set of Feynman rules for Yang-Mills amplitudes based 
on scalar propagators connecting three and four point vertices. The starting point 
is the expansion of the axial gauge propagator in terms of polarisation vectors, 



' 'et{k)e-{k) + e-^{k)et{k) + el{k)e1{k)], 



(3.1) 



k 

where 

e^(k)-^^^ e-(k)-^^hA!ll e^(k)-2^q (3 2) 

Here g is a null reference momentum which may be complex. For any momentum k 
we define its q-nullified form 

Contracting the polarisation vectors into the usual Yang-Mills three-point vertex 
yields the familiar three-point MHV and MHV vertices, 



iv^ " • ' (23){31) |1,||2,] 



tV2 " • • ' |23||31| {1,)(2«)' 

along with a V3(l^, 2~, 3*^) vertex. In the formula above, all momenta are g-nuUified. 
As vertices of this last type must be attached together in pairs, it is natural to 
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absorb the resulting four-point configurations into an effective four-point vertex along 
with the Yang-Mills four-point vertex. These effective four-point vertices contain 
prefactors 

[pq] (mg) 



and 



(3.5) 



(pq) [mq] ' 

for each positive-helicity leg p, and each negative-helicity leg m, respectively. 

When adopting a recursive approach which involves shifting a negative-helicity 
leg a and a positive-helicity leg b, the recursion-optimised choice for the reference 
momentum q is 

Xg = Xa, Xg = Xf,. (3.6) 

With this choice of q the prefactors of four-point vertices ( p.5|) involving a shifted 
leg vanish. Furthermore from (p.4|) the legs a and b can only enter a diagram on an 



MHV or MHVthree-point vertex respectively. 

Thus for the single-minus amplitudes, at tree and one-loop level, the external 
negative-helicity leg must enter the diagram via an MHV three-point vertex and this 
must have a negative-helicity internal leg. This leaves insufficient negative helicities 
to have a four-point vertex anywhere in the diagram. At tree level there are no 
non-vanishing diagrams whilst at one-loop we have a single MHV three-point vertex 
and several MHV three-point vertices in each diagram. These rules apply to both 
Yang-Mills and gravity calculations. For gravity we define the tree amplitudes using 
the Kawai-Lewellen-Tye (KLT) expressions [36]. 





Figure 3: With the constraints that (1) the negative-helicity leg enters via an MHV 
three-point vertex and (2) the four-point vertices vanish, we only have non-vanishing 
diagrams with a single three-point MHV vertex with the remaining vertices three- 
point MHV with internal helicities organised as shown in these sample diagrams. 



We now wish to characterise the singularities when either Sbc or Say vanish. 
Singularities arise in the integration from the region of loop momentum where the 
denominators of three adjacent propagators vanish simultaneously, as the two null 
legs to which they connect become collinear. The diagrams of interest for any single- 
minus amplitude can then be collected into the forms shown in figure ^ Note that 
we evaluate these diagrams for real momenta and only carry out analytic shifts on 
the final expressions. The circles in these diagrams represent the sums of all possible 
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tree diagrams with two internal legs and the given external legs. We denote these by 
r(a, b, . . .). In the integration region of interest all the legs of r are close to null and 
the internal legs are close to coUinear. 




Figure 4: Singularities in Sbc and Say arise in integrations over the terms shown. 



In each case there are three options for the helicities within the loop, as illustrated 
in figure ^ Let us consider figure With the configuration of figure r vanishes 
at the integration singularity because it is a one-minus tree amplitude in this limit 
and so the diagram has vanishing residue. 




(a) 





(c) 



Figure 5: The three possible helicity structures of figure 4a. 



The diagram 5b evaluates to 



[b\l\a)[c\l\a) {l-c, ayT{{l-c)+,d+,--- ,a-,(/ + 6)-) 



{ba){ca) {l + b,ay 



Pii + hYii-KY 



(3.7) 



where the momenta in the spinor products are g-nullified as in (|3.3|) . We construct 
a basis for the loop momentum built on b and c via 

I = ai{h + K) + Oi2{h - kc) + ("3 + ioii)'^jz~\^bK + ("3 - ioLi)T—T^c\ (3.8) 

{ba) {ca) 

Under this parametrisation, 



fm = 1. 



dai F{ai)f{l{ai)) 



(3.9) 
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where F{ai) has no dependence on Sbc ^■ 
Also, 

[b\l\a) = (ai — 0:2 + «3 + iai) [b\c\a), 
[c\l\a) = («! + a2 + as — 10:4) [c\b\a). 

After these manipulations the integrand from figure ffi becomes 



(3.10) 



Wi+f^^^^^ - c)\ d\...,a~,{l + b)-) X F'{a,). (3.11) 

When /, b and c become coUinear, r approaches the coUinear limit of an MHV tree 
amplitude. Within r there are diagrams with an explicit Sbc pole. The singular factor 
from the integration around the coUinear limit combines with the explicit pole factor 
to give the double pole discussed previously. In addition to the leading behaviour of 
r in the coUinear limit, we need to know its finite piece in order to determine the 
residue at this pole. The diagram |5^ gives the same contribution. 

We can apply a similar analysis to the contributions from diagrams of the type 
shown in figure In this case r approaches either a one-minus or an all-plus tree 
amplitude in the region of interest and so vanishes. Thus diagrams of the type shown 
in figure ^ give no contribution. 

In order to evaluate the contribution from ( p.ll|) we must evaluate the tree struc- 



tures to order (bc)^. For diagrams within r involving l/s^c this means going beyond 
leading order. The loop part of these diagrams is a triangle and the calculation is 
readily done exactly. The diagrams without this propagator need only be calculated 
to leading order. In this regard, not only is the recursive approach selecting a subset 
of diagrams for calculation, it is also allowing us to calculate these diagrams in a 
very convenient limit. 

In the following section we apply an augmented recursive analysis to the calcula- 
tion of the amplitude M1-'°°p(1-, 2+, 3+, 4+, 5+). For gravity, using the KLT relations 
for tree amplitudes [36] the equivalent expression to ( ^.11| ) is 



^%av((/ - c)+, d+, . . . , a-, (/ + b)-) X F'(a,). (3.12) 
{oc) (/ + b, a) 



4. The graviton scattering amplitude M1-i°°p(1 ,2+,3+,4+,5+) 

To compute this amplitude recursively, as discussed in the previous section, we must 

^We have not explicitly introduced a regulator, although even for the finite amplitudes under 
consideration individual diagrams may diverge. If we were to use a Pauli-Villars regulator with 
mass Mpv for instance, we could still extract the same momentum-dependent prefactor but the 
remaining integral would depend on the and Mpy/shc- Knowing that any divergent pieces cancel 
in the full amplitude allows us to consider only the finite pieces, which are independent of Afpy- 
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compute three types of contribution: 





l-c , 


/ d+ 


@" 


- e+ 


l+b^ 





(1) (2) (3) 

together with the contributions obtained by summing over the distinct permutations 
of c, d and e. 

The first two of these involve single poles only, so we only need the loop structures 
to leading order and we can use the corresponding four-pt one-loop amplitudes. The 
final structure contains a double pole so we must evaluate both the tree structure on 
the right and the loop pieces more carefully. The first diagram uses the four-point 
one-minus amphtude whereas the second requires the four-point all-plus amplitude, 
both are given in ( [B.4| ). We obtain, 

1 (ad)^(ae)^[6c][rfe]^ {{cd)^ae)^ + (ac) (cd) (de) (ae) + {ac)^de)^) 
R,[a,b,c,d,e} - ^ (a6)2(6c)(ce)2(crf)2(rfe)2 

i?2(o,6,c,d,e) = -^pj^^l^ {[bcnde]' + [bc][cd][de][be] + [cdnbe]') . 

(4.1) 

For the third diagram, we need the tree diagrams which constitute Tgrav of ( |3.12|) . 
Mindful of the recursive analysis that we will ultimately perform, we calculate these 
diagrams as Laurent series in (be), dropping terms that will not contribute to the 
residues. 

We require the five-point contributions with two off-shell legs B~ and C^, car- 
rying momenta B = 1 + b and C = c — /, respectively. 

The KLT relation [36] between Yang-Mills amplitudes and gravity amplitudes 
at five points is 

M{a-, C+, d+, e+) =SBcSdeA{a- , B' , C+, c/+, e+)A(a", C+, B' , e+, 

+SBdSceA{a', B-, ci+, C+, e+)A(a-, d+, 5", e+, C+), ^^'^^ 

where we have chosen a form of the KLT relations that restricts the (be) pole to the 
first term. The KLT relations are only valid for on-shell momenta, although these 
momenta may be in higher dimensions. If we assume the deviation from eq. (|4.2| ) 
may be neglected in the region around 5^ = = 0, we see that the gravity tree 
structure has the form, 

(T^lcading \ /T^leading \ 
i_ T^sub-leading \ I , T^sub-lcading \ /* o\ 
{be) ^ J\ {be) J' 
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where all diagrams contribute to the sub-leading pieces but only diagrams involving 
a Vs{B~ , C"*", x) vertex contribute to the leading pieces. The second term in eq. (|4.2| ) 
is only needed to leading order and its contribution to the residue will be directly 
determined by the on-shell Yang-Mills MHV amplitudes. The amplitude generated 
using the leading and sub-leading singularity terms from ( ^.2| ) has the correct sym- 
metries and coUinear limits. Additionally, the five-point amplitude has been verified 
by a completely independent string-based rules computation. The general case is 
worthy of further study [37]. 

First we establish the double pole term. This arises from the poles in each of the 
Yang-Mills tree amplitudes in the first term of ([4.2|) . We evaluate this diagrammat- 
ically. The Yang-Mills amplitude, A{a~, B~, C~^, d~^, e^) receives contributions from 
five diagrams. The two which contribute to the pole are: 



J- h- 



B- a- B- 

(a) (b) 



a 




with contributions 



{Baf {a\hc\a)[de][eh] 
{Ca)"^ Sbc{da) (ea) [ae] [ab] 
(Ba)'^ {a\bc\a){ca) 



{CaY {bc)[ab]{de) {da) {ea) 



fa{ai), 



(4.4) 



where we have used {a\BC\a) = {a\{b + l){c — l)\a) = {a\bc\a) f {ai) , etc. The param- 
eters contained in and are the same for both diagrams and the sum of the two 
contributions is 



Da + Db 



{Ba)^ {a\bc\a) ( [b\ad\e] - [b\cb\e\ 
{Ca)"^ Sbc[ob]{da) {ea) 



[ae] {de) 



'fa{ai), 



(4.5) 



where the second term is sub-leading in the {be) pole. 

The leading pole in the other Yang-Mills factor is obtained analogously and, 
combining, we obtain the leading pole in (|4.2|) , 

{Ba)^ 



^Sbc^de' 



{ab) {ac) [de\ {ab) {ac) [de\ 



-f>^)■ 



(4.6) 



{Ca)^ {be) {ea)[ae]{de) {be) {da) [ad] {de)' 

Combining this with the factors arising from the left hand part of the full diagram 
and integrating over the the leading term in the Laurent series is proportional to 

[6c] 3 



{be) 



[be][de] 



{ab){ae)[de] {ab){ae)[de] 
{ea) [ae] {de) {be) {da) [ad] 



(4.7) 
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which clearly displays the double pole factor. The constant of proportionality is most 
readily fixed by looking at coUinear limits. 

We must now enumerate the contributions that are sub-leading in the s^c pole. 
These come from a variety of sources. We express these single-pole terms as the 
double-pole factor of (^T?!) , multiplied by a factor 6. Firstly, we have the sub- leading 
contribution of ( |4.5|) together with the corresponding contribution from the other 
Yang-Mills factor, 

^ Sbcjbe] Sbcjbd] 
' [b\ad\e] [b\ae\dy ^ ' 

Next we have the sub-leading diagrams for the Yang-Mills amplitudes in the first 
term of ( [4.2|) shown below: 

B- Ct d+ B- 

ki-+ ks 






a" e+ a" ^ 

(c) (d) (e) 

We note that in the first two of these the ki propagator feeds into the two diagrams 
that would make a one-minus four-point tree if ki and C were both null. As we know 
that this vanishes when all the legs are null, the sum of the first two diagrams must 
be of the form: C^X + kfY. We can drop the terms containing a factor as we are 
already at sub- leading order, leaving something proportional kf. Thus taking both 
terms together leads to the cancellation of the SaB propagator: 

^ ^ (Ba)"^ [b\B\a) (ca) „ . , ,^ ^, 

{Ca)^ [ab\{ea) {de) [cd) 



(4.10) 



Pulling out a factor of ( |4.7|) leaves 

^ _ Sbc[e\a\c) 
^ Sab[e\d\c)' 

We can apply the same procedure to the final diagram giving 

{bc){de) ne\B\a)[eb] ^ [d\B\a)[db] \ 
■Sab[c^e] \ {da){cd) {ea){ce) J 

Finally we need the second term in (|4.2|), SBdSceA{a, B, d, C, e)A{a, d, B, e, C), which 
we evaluate using MHV tree amplitudes. After extracting the double-pole factor we 
obtain 

_ {bc){de)[d\B\a)[e\C\a) 
^ [bc][de]{aby{cd){ce) ' ^ ' ' 
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We thus have the leading and sub-leading poles expressed as 

H[6c]Me]MM?JM^^ xfl + Ts] (4.13) 

(6c) ^ ' {ea)[ae]{de) {bc){da)[ad] \ ^ ) 

We can now use the pole to determine the amplitude recursively. This involves 
applying the shift ( |1.1| ) and evaluating aX z = —{bc)/{ac). The coefficient of the 
double pole in ( [4.7| ) has a z dependence under this shift which generates a further 
contribution to the single pole since 



Res — 

,2(Z — Zj)^ / Zf Zi dz 

Carrying this out and combining with the contributions of the 5s gives 

I {ad) {he) l (ae)(6c) 
^ ' ' ' ' ^ 2(a6)(crf) 2{ah){ce) 

^ {dh][eh] {bc){de) [dc][ec] {be) {de) {ca)' 
{dc){ec) [be] [de] {dc){ec) [be] [de] {ba)^ ^ ' ' 

_ 7 [d(^[eb] {be) {de) {ca) _ 7 [db][ee] {be) {de) {ea) 
~ 2 {de) (ec) [be] [de] {ba) ~ 2 {de) (ec) [be] [de] {ba) ' 

The full one-minus amplitude can now be written as the sum over recursive 
contributions arising from three orderings of the external legs, 

M^-'°°P(1-, 2+, 3+, 4+, 5+) = R{1, 2, 3, 4, 5) + R{1, 2, 4, 5, 3) + R{1, 2, 5, 3, 4). (4.16) 

with the full amplitude having a factor of zk^/IOtt^ as in ( p.l|) . 

Each recursive term is a sum over the three classes of recursive diagram, 

R{a, b, c, d, e) = Ri{a, b, c, d, e) + i?2(a, b, c, d, e) + Rs^a, b, c, d, e), (4-17) 

where Ri and R2 are given by ([4.1|), and 



„ / , , N 1 {ab)^ {ae)'^[be]'^[de] , . . , , x\ /. 

6, c, d, e) = + A(a. c. d. e)) . (4.18) 

The overall normalisation can be obtained by evaluating the parameter integrals or, 
more easily, fixed by f actor ising the known four-point amplitude. The individual 
factors on the terms in A are also obtainable by parameter integration or more 
conveniently by the normalisation of the coUinear limits. 

This form for the amplitude has the correct collinear limits and is symmetric 
under interchange of any pair of positive-helicity legs. We have also checked that the 
amplitude agrees with that calculated by string-based rules. This calculation can 
readily be extended to the six-point case, M-^"^°°p(1~, 2+, 3"*", 4+, 5+, 6"*"). We have 
constructed the amplitude and again checked that it has the correct symmetries 
and collinear limits. This result is presented in appendix 0. Mathematica code for 
both the five- and six-point amplitudes may be found at http://pyweb.swan.ac. 
uk/~dunbar/ graviton . html. 
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5. Conclusions and remarks 



In this article we have demonstrated how to augment recursion in order to determine 
the rational terms in amplitudes with double poles under a complex shift. Double 
poles are generic in amplitudes, however it is often possible to carry out a recursion 
which avoids then. However, double poles are unavoidable in the case of the one- 
minus Yang-Mills amplitudes y4^"'°°P(l~, 2+, S"*", . . . , and the gravity amplitudes 
M^"'°°P(1~, 2+, 3^, . . . , n"^). In the absence of a universal soft factor analogous to 
( |1.8|) , in order to perform the augmented recursion the sub-leading poles must be 
determined on a case-by-case basis. While we have done this for both the five- and 
six-point one-minus gravity amplitudes, this procedure could be used to calculate 
the seven-point or indeed any higher-point one-minus amplitude. 



A. Six-point single-minus amplitude 

The six-point one-loop single-minus graviton scattering amplitude can also be cal- 
culated using augmented recursion. The calculation follows that of the five-point 
amplitude with the addition of factorisations involving a four-point tree amplitude 
and a four-point loop amplitude. The shift employed is once again Ai — Ai — zX2^ 
A2 — )■ A2 + z\i. The amplitude is given by 

M^-'°°P(1-, 2+, 3+, 4+, 5+, 6+) = 

xg{3,4,5,6} d(6)/o h ^\ 

{s/l,y2,»/3}U{a;}={3,4,5,6} + -"-3 l^; 3;| i, ?/2, 2/3i ^a; j / 

+ i?f^(l,xi,a;2|2,?/i,?/2; 2;xi,a;2). (A.l) 

{i'l,a;2}C{3,4,5,6} 
{yi,y2}U{a;i,X2}={3,4,5,6} 

In each of these terms the vertical bar denotes a split of the external momenta with 
the relevant pole arising when the shifted total momentum to the right of bar is null. 

The terms are the factorisations involving a three-point MHV tree and a 
five-point all-plus one-loop amphtude: 

i?f)(l,x|2,yi,i/2,z/3;^.) = % I! 1 M^-'°°P(p+, 2+, y+, y+, ^3+), (A.2) 

[12] [Ix] 

with = [1 x] / [2 x] and p = (A^ + A^ [1 2] / [x 2])A^. Similarly, the R^^ terms are the 
factorisations involving a three-point MHV tree and a five-point one-minus one-loop 
amplitude: 

Rf\2,x\l,y^,y2,ys;z.) = S |! M^-'°°p(1-, y+, y+, ^3+, (A.3) 
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with = —{x2) I and p = A^'(A^ + (12)A^/ (l^c)). The J^^^ terms are the 
factorisations involving a four-point MHV tree and a four-point all-plus one-loop 
amplitudes: 

i?^ ^(1, xi, X2I2, yi, y2; ^0:1x2) = 7 > (A.4) 

''1X1X2 



with z^;^^;^ = iixix2/[2|-Pxix2|l) and p = /Cj^^ /Ca^^- 

The i?3^^ terms are the augmented pieces arising from the (2x) poles. Here, 
z,^-{x2) I {xl). 

Ylx^ { 

R3 (2,a;|l,7/i,7/2,|/3;^a;) = | -KLTpll, 2, |/2, l/s) 

+ -^S3/2j/3YMl(1, 2, 1/1, 1/2, 1/3) 

X Kri/3YML(i,lx,y3,yi,y2) + (sj,ij/2 + Sj/i2/3)YML(i, 2, x, ^3, y2, Z/i)] 
+ [22;]sy22/3YML(i, 2, yi, ^2, Vs) 

X [sj,,j/3YMs(i, 2, x, ys, yi, ^2) + (sj/12/2 + Sy,ys)YMs{i, 2, ^3, y2, yi)] 
+ [2a:^]s2/22/3YMs(i, 2, x, 7/1, ^2, ^3) 

X [sy,y^YMi^{i, 2, y3, yi, Z/2) + + Sj/i2/3)YML(i, 2, ^3, ^2, yi)] (A.5) 

\2x] - - 

+ ^2^Sy^y^YMj^{l, 2, ^/2, yi, 7/3) 

X K2j/3YML(i,2,x,y3,y2,yi) + +%22/3)YML(i,2,x,y3,yi,y2)] 
+ [22;]s2/i2/3YML(i, 2, 7/2, yi, ys) 

X [s2/2j/3YMs(i, 2, X, ys, y2, yi) + (sj/22/1 + S2/22/3)YMs(i, 2, ys, yi, ^2)] 
+ [2a;]s2/i2/3YMs(i, 2, ^2, yi, y3) 

X [*2/2j/3Y^l(-'-' 2, a;, 1/3, 1/2, yi) + {sy^y^ + s^j^g )YML(l,2,a;,|/3,|/i,|/2)] 
where the KLT terms contributing to the double pole have leading Yang- Mills factors: 



(n h n ^ ^ f\- (ac)(&a) \ [hc]{ca) [ [f\Pde\a) , [d\Pef\a) 
' ' ' ' H(rfa)(ea)(/a)\ U^f \ (de) + (/e) 

[d|Pe/|a)[&|Pe/|«) , [d\Pef\a)[ef][fb] [f\PM[fb] \ 
{ef)Uf, ^ [af]tefa [af]{de) J' ^ 
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and sub-leading factors: 

(ha) J [./|P,;|«)[6|P,;|a)^ 



YMs(a,6, c,d,ej) 



[ab]{da){ea){fa)\ {ef)tlf^ 

^ [d\Pef\a)[ef][fb] r mf\a) ^ [bfV\ _ [me\a)[bff 

[af]tefa V Wa [«/] / [af]'^{de) 
^ 1 [bc]{acy _ l [f\h-h\a){ac)[fb] 

2[ab]{cd){de){ef){fa) 2 [ab]{cd){de){ea)sfa 
_ l {a\ih-K)Pef\a){ac) nb\Pef\a) ^ [bf][ef] \ 

2 [ab]{cd) (da) {ea)tefa \{fe){fa) Sfa J' 



(A.7) 



Finally the finite terms in the KLT sum are: 

KLTF(a, b, c, d, e, /) = 



s. 



,f{ab) 



4 



2{ad){ae){af){bdy{be)^bf){ef)^ 

X \ {ae){6{a\bPef\b)[d\b\a) + 7{a\bPef\b)[d\c\a) + 7{a\cPef\b)[d\b\a) + 6{a\cPef\b)[d\c\a)) 



+ 



+ {ab){6{e\fb\a)[d\b\a) + 7{e\fb\a)[d\c\a) + 7{e\fc\a)[d\b\a) + 6{e\fc\a)[d\c\a)) 
Sdf{aby 



2 (ad) (ae) (bd) (be) {bf)^ {de) {df) 

X {6[e|6|a)[/|6|a)+7[e|6|a)[/|c|a) + 7[e|c|a)[/|6|a) + 6[e|c|a)[/|c|a)} 



+ 2MRMMMMW^^^^ ^ M|%)[e|c|a)[/|%)} 

{d O e} 

(A.8) 

Expressed naively, without attempting optimisation, as a rational polynomial 
of the it has a LeafCount of 355,053. For comparison, the LeafCount of the 
five-point one-minus gravity amplitude is 4,549, and for the six-point one-minus 
Yang-Mills amplitude is 1,541. 

B. Graviton scattering amplitudes 

We define tree and one-loop amplitudes in gravity for which all field couplings have 
been removed, i.e., 

M'ril, 2, . . . , n) = i^^^'^-^^Mr (1, 2, . . . , n), 
>l^-^-P(l, 2, . . . , n) = ^M^-i-P(l, 2, . . . , n). 
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As for Yang-Mills amplitudes we express amplitudes using the spinor helicity for- 
malism. For the four dimensional case there are only two graviton helicities and 
their polarisation tensors can be constructed from direct products of Yang-Mills 
polarisations vectors, 

= <4, £^. = £^£u- (B.2) 
If we consider the Feynman diagrams for a gravity one-loop scattering amplitude, 
performing a Passarino-Veltman reduction [38] allows us to reduce any one-loop 
amplitude to the form 

M^-i°°P(l, ...,n) = J2cJl + Y^d,Ii + J2ekI^ + R+ 0(e). (B.3) 

i j k 

Relatively few graviton scattering amplitudes have been computed. In fact, only 
the four-point amplitudes have been computed for all helicity configurations [31,32, 
39-41] and all possible matter types circulating in the loop. For four points there arc 
three independent helicity configurations for the external gravitons: M(l+, 2+, 3+, 4+), 
M(l~, 2+, 3^, 4+) and M(l^, 2^, 3"'", 4"*"). The all-plus and one-minus vanish at tree 
level and have one-loop amplitudes which are purely rational (to order e°). These 
amplitudes, for any matter content, are 

^ ' '"^ ^'v(12)(23)(34)(41); 1920 ' 

^ ' ' ' ^ '\u) V[12](23)(34)[41]/ 5760 ' 

where s — {ki + k2)'^, t = (A;i + /c4)^, u = {ki + k^y and Ng = Ns — Np is the number of 
bosonic states in the loop minus the number of fermionic states. The amplitudes for 
pure gravity are found by putting A^^ = 2 in the above expressions since a graviton 
has two helicity states. These amplitudes vanish in any supersymmetric theory. The 
n-point all-plus and one-minus amplitudes are also particle-type independent up to 
a prefactor of A^^, as can be seen from the vanishing of these amplitudes in any 
supersymmetric theory as a consequence of supersymmetric Ward identities [42,43]. 
It is therefore sufficient, for these configurations, to compute the amplitude with a 
scalar particle circulating in the loop. 

Beyond four points most of the explicit graviton amplitudes are for scattering in 
supersymmetric theories. For J\f — S supergravity the n-point MHV is known [44] , as 
are the NMHV six- [45,46] and seven-point amphtudes [47]. In ref. [44] a 'dimension 
shift' relation [24] allowed the conjecture of an ansatz for the all-plus n-point ampli- 
tudes. This amplitude is an ingredient in the recursion of the one-minus amplitude. 

C. string-based rules calculation of M1-i°°p(1-, 2+, 3+, 4+, 5+) 

The string-based rules were introduced in refs. [29,30,48] as a method of calculating 
(one-loop) gauge theory amplitudes. Their extension to gravity, in the form we 



-17- 



use here, was given in refs. [31,32]. In this appendix we summarise these rules 
and then describe how they are apphed to compute M^"^°°p(1~, 2"*", 3"*", 4"*", S"*"). Our 
presentation treats the method as something of a 'black box' for obtaining field- 
theory results and we refer those interested in the details of its string-theory origins 
and derivation to the literature. 

C.l Summary of the string-based rules for gravity amplitudes 

String-based rules use one-loop (p^-like graphs to compute the one-loop corrections 
to a field theory amplitude. The terms produced take the form of a rational func- 
tion of the kinematic variables, within a Feynman parametrisation of a tensor loop 
integral. This approach has the advantage of significant computational savings over 
the traditional Feynman graph method: far fewer graphs are involved and early ap- 
plication of simplifications from the spinor-helicity formalism reduce the complexity 
of the associated expressions. 

We begin by drawing all one-loop 0^ graphs excluding massless bubbles (which 
vanish in dimensional regularisation) and tadpoles. We label the outermost legs of 
the graphs with the particles' momenta, /ci, . . . , An internal line bears the same 
label as the first line or leg found going anti-clockwise about its outermost vertex. 
(For examples of such graphs and labellings, see figure ^.) All independent labellings 
of external legs contribute for gravity amplitudes. The one-loop correction to the 
amplitude is then given by 

Mi-i-P(l,2,...,n)= ^(^)' 

graphs 7 

where the contribution from a graph 7 with an n^-propagator loop is 



{'^l<k<l<nt ' -Pii^ikili^ ^ikil)} 

(C.l) 

In this formula, ii, . . . ,in^ are the labels of the lines adjoining the loop going clock- 
wise. Xij = Xi — Xj, with Xi„^ fixed at 1, and Pj^, is the momentum entering the 
loop along the line with innermost label ik- The 'reduced kinematic factor' for 7, 
K^"'\xij^, . . . ,a;j^^_J, is a polynomial in the Xi^ which for a gravity theory with no 
supersymmetries is of order 2ne. 

We compute K^^'^Xij^, . . . , Xi^^_^) as follows: the starting point is the overall 
graviton kinematic factor 



Y\ dxidxi JJ^ exp jfcj ■ kj G% + {hi ■ ej - kj ■ ei)G'^ - Ei ■ ej 

i=l l<i<j<n 

+{ki-ej-kj-ei)G'i-ei-ejd'i-{e,-ej-ej-ei)H'^\ , (C.2) 

J multi-linear 
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where 'multi-linear' indicates that we retain only the coefficient of YYi=i ^i^i 
expansion of the exponentials. The graviton polarisation tensor are then recon- 
structed using ( [B.2| ). K. contains much structure from the string theory perspec- 
tive: = GB{xij) is the bosonic Green's function on the string world-sheet, and 
the Xi{xi) are closed-string left(right)-moving co-ordinates. Other objects present 
are the derivatives of G^: Gg = dG^^/dxi, which is antisymmetric in i,j, and 
Gg = d'^G^^/dxf (with similar expressions for the right-moving Gg and Gg); and 

= d'^G^^/dxidxi. However, for our purposes we will simply treat (|C.2|) as an ob- 
ject for obtaining reduced kinematic factors by the application of some substitution 
rules that implement the field theory limit. 

For the helicity configuration under consideration and a judicious choice of ref- 
erence momenta for the polarisation vectors, we shall see that the coefficients of the 
second-order derivatives of the Green's functions vanish. Nevertheless, in general 
this is not so and we should eliminate the Gb and Gb from /C using integration by 
parts, which may lead to additional factors of the appearing. Each factor 
should then be eliminated by replacing it with the Feynman denominator relevant to 
the diagram under consideration [i.e. the expression found within the curly braces 
in ( |G.1| )). At this point we simply drop the integration over the world-sheet co- 
ordinates and the leading factors of exp(A;j ■ fej-Gg) from K (their contributions are 
built into the rules). 

Now consider a consecutive pair of lines joining on to the loop in a (j)^ graph, 
labelled («, j) going clockwise. We can 'pinch off' this pair of lines by attaching them 
instead to a new vertex and then drawing a new line carrying the label j from this 
vertex back to the loop. Any of the graphs drawn for string-based rules may be 
obtained this way (for example, the graph of figure |6^ is obtained from figure |6a|ffist 
by pinching off (2,3), and then (3,4)). For each such pinch used to reach a graph 
we: (1) discard all terms in the expression obtained above except those containing 
exactly one power of GgGg; (2) replace i with j in all remaining Gb, Gb; and (3) 
multiply by —1/kfj, where kij is the momentum carried by the new line formed by 
the pinch. 

Next we apply the substitution rules. These act on the derivatives of the Green's 
functions in a left /right-independent manner, replacing them with polynomials in 
the Xij in a way that depends on the particle content of the loop. In particular, 
the rule for a single scalar degree of freedom running around the loop is the simple 
substitution 

Gb, Gb Xij - - signxij. (C.3) 

There are other rules for particles of higher spin in the loop (including the graviton), 
but by the discussion in appendix |B|, (|C.3| ) is all we need for a one-minus amplitude. 
Therefore we compute the amplitude by applying ( |G.3|) to the reduced kinematic 
factors and multiplying by A^^ = 2. Finally we make change of the integration 
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variables in ( |C.1| ) to ak using Xi^, = Yl'i=i^i^ which yields an integral in the usual 
Feynman parametrisation. 

C.2 Application to M1-'°°p(1 ,2+,3+,4+,5+) 

The (topologies of the) graphs that have a non- vanishing contribution to M^"^°°p(1^, 
2+, 3+, 4+, 5^) are shown in figure There are 117 such labelled graphs in total: 12 
massless pentagons (figure |6a|) , 30 one-mass boxes (figure 15 two-mass triangles 
(figure 0) and 30 one-mass triangles (figure There are also 30 massive bubbles, 



but these vanish by the pinching process when using the spinor helicity choice ( |C.4| ) 
below. 





(d) 

Figure 6: Topologies for 0^-like Feynman diagrams that have a non-vanishing con- 
tribution to the string-based rules calculation of M^"^°°p(1~, 2"*", 3^, 4+, 5+). (The 
labellings shown are non-vanishing examples; other orderings also contribute.) 



In order to define the Ei and Si (which are set to the same values after multi- 
linearisation in ( p.2| )), we choose fcs as the reference momentum for the first graviton 
and ki for the rest. In the spinor- helicity formalism the polarisation vectors are 



[5|7ll) 
V2[l 5] 



and 



for i ^ 1. 



(C.4) 



We have the standard spinor-helicity results that ki ■ Si = ■ ei = ki ■ Si = for all 
i, and furthermore for this choice Ei ■ Ej vanishes for all i,j, so there are no second 
derivatives of Green's functions to handle. After dropping the exp{ki ■ /c^-Gg), (^ 
becomes, 
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X {k2 ■ e,G'^ + k, ■ esGl^ + h ■ e^G''i){k2 ■ e.G^i + k, ■ e.G^i + h ■ e.G^i) 

X {k2 ■ £5 + fc3 ■ £5 G^i + k,-e, x (1 ^ r). (C.5) 

Here, '(1 — )■ r)' denotes taking the expression to the left and replacing all Gg with 
Gg. We can now see why there are no bubble graphs in the problem. They come 
from three pinches that form two independent trees, but any such sequence of pinches 
will either pull out a factor of the form {ki + ki) ■ Ei, which vanishes by conservation 
of momentum and the remarks below ( |C.4|) , or simply run out of pinchable Gbs. 
Since the same substitution rule (|C.3| ) is applied independently to both the left- 



and right-moving sectors, and at each step in the pinching we pull out terms contain- 
ing exactly one power of both Gg and Gg, we can in fact proceed in a rather more 
straightforward manner by applying the pinching and substituting for just the left- 
moving factors of ( |U.5| ), then taking the square of the result as K^'^\xij^, . . . , J, 
taking care not to square the kinematic factors that arise from the trees during 
pinching. 

We do this for all 117 graphs and substitute back into (|C.1D , changing the vari- 
ables to the usual Feynman parameters. Each graph thereby yields an expression of 
the form 

P(7) = J]xW(ri, . . . , r„J/W[a? ■ ■ ■ aZ% (C.6) 

fe} 

where X^^'^ri, . . . , r„J is a rational coefficient. The n^-gonal tensor Feynman integral 
with momentum configuration relevant to the graph 7 is defined as 

/£^[a?---n7]=rK-2 + e) fd-'a " '"j^^^ ' ^l^l,,^ , (C.7) 

with the array Sl]^ given in terms of the momenta Pi entering 7's loop by 



' ioT k = l, 

i(Pj^ H hPi(„J^ otherwise. 



0(7) 



These integrals may be evaluated by the recursive approach detailed in refs. [34, 
49]. In this way we have constructed an integral database using computer algebra, 
indexed by the tuples {(ri, . . . ,r„J| ^j^^i < 2n£} for integrals with 3 < < 5 and 
up to 5 — rii massive legs. The table contains both the rational pieces of the integrals 
and the rational coefficients of their (di) logarithms. Since M^"^°°p(1~, 2~^, 3~^, 4+, 5^) 
is entirely rational, we can use the vanishing of the logarithmic terms as a consistency 
check. Computationally the most complicated integrals are the pentagons where we 
must evaluate integrals with Feynman parameter polynomials of order ten. 
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The coefficients produced in this manner are far too large and cumbersome 
to present here (or even form a compact analytic expression for the amplitude at 
present); nevertheless, they are amenable to exact numeric arithmetic at a kinematic 
point. The results have been compared with, and agree with, the recursion-derived 
expression ( ^.IGp . 
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